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Abstract. This brief review describes radio observations of colliding winds 
in massive stars starting with the first direct observational support for the 
colliding-wind model advanced in the early 1990's to explain non-thermal radio 
and thermal X-ray emission in some massive stars. Studies of the well-studied 
and highly-eccentric WR-l-0 star system WR 140 are described along with recent 
observations of 0-star systems. Also discussed is the binary nature of almost all 
massive stars that exhibit non-thermal behavior and some strategies for finding 
new systems. 



1 Introduction 

Massive stars have dense stellar winds that are photo-ionized by the strong UV 
radiation fields from the underlying star, giving rise to an excess of free-free 
continuum emission over the photospheric emission, that can be observed from 
IR to radio wavelengths. The brightness temperature of this emission is ~ 10^ K, 
as expected from a photo-ionized envelope. The emission is partially optically 
thick, and between mid-IR and radio frequencies has a power-law spectrum Su oc 
z^", with spectral index (a) typically ~ -1-0.7 (e.g. Williams 1996). A number of 
massive stars have radio emission properties that differ from this typical picture: 
they exhibit brightness temperatures ~ 10^ — 10'' K, and flat or negative spectral 
indices, both properties characteristic of non-thermal synchrotron emission and 
therefore of relativi s tic el ectrons. Often, this non-thermal emission is variable 
(|Abbott et al.lll984l . Il986l 1. 



It was first noted by Ivan der Hucht et al.l (|l992l ^ that many Wolf-Rayet 



(WR) stars exhibiting non-thermal emission were long-period binary systems, 
and they proposed that the non-thermal emission arose from the region where the 
massive stellar winds of the two binary components collided, the wind-collision 
region (WCR) outside the optically thick part of the circumbinary stellar wind 
envelope . The colliding-wind b inary (CWB) model was described in the seminal 
paper of lEichler fc Usovl (|l993l 1. 



Spatially resolved observations of the WR-l-0 system WR 147 (jWilliams et al. 



19971 ) presented the first unequivocal confirmation of the CWB model. The 



model is supported further by the dramatic variations of the synchrotron radio 
emission in the 7.9-year WR-l-0 binary system WR 140, that are clearly modu- 
lated by the binary orbit (c.f. Fig. [3]). Among WR stars, the CWB model has 
now been substantiated with over 90% of WR stars that exhibit non-thermal 
emission are either a bina ry system or have a "nearby" massive stellar compan- 
ion that results in a WCR ( Dougherty &: WilliamslfioOOl ) . Amongst 0-type stars. 



there is now increasing evidence that a similar explanation also holds (see papers 
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by van Loo, and Benaglia in these proceedings; De Beckei] ( 20071 )). The nature 
of the acceleration mechanism is, as yet, unidentified but is widely attributed to 
diffuse shock acceleration in shocks that "bound" the WCR. Alter native mecha- 
nisms, such have magnetic reconnection have been discussed (e.g. iJardine et al 
19961 1 



Colliding-wind binary systems present an important laboratory for investi- 
gating the underlying physics of particle acceleration, of not only the relativistic 
electrons that give rise to the observed radio emission, but potentially to ion ac- 
celeration that may give rise to gamma-ray GeV-TeV emission (see other papers 
in these proceedings) . In addition to the excellent modelling constraints provided 
from high-resolution VLBI imaging, radiometry and most recently high-energy 
observations, CWBs provide access to higher mass, radiation and magnetic field 
densities than found in supernova remnants, that have been widely and success- 
fully used for studying particle acceleration, particularly diffusive shock acceler- 
ation (see Ellison, these proceedings). 

Here the observational evidence for non-thermal radio emission arising in 
WCRs in both WR and 0-star binary systems is reviewed. In closing, some 
potential schemes for finding more of these particle acceleration laboratories are 
suggested. 



2 Observational evidence for colliding-wind systems 

2.1 The first essential proof of a CWB - imaging WR147 



Williams et al. ( 19971 ) provided the first direct imaging confirmation of the collidin g- 



wind model in observations of WR147. This system was shown bv lMoran et al 



(119891 ') using the MERLIN telescope array in the UK to have two distinct com- 



ponents; a thermal component associated with the optical positio n of the WR 
star, a nd a non-thermal component 0.6 arcseconds to the North. iMoran et al 



( 19891 ) speculated that this may be due to the presence of a companion star 
as had been suggeste d previously to explain non-thermal emission in WR140 
(IWilhams et al.l ll987l 1. No companion could be found in a search for objects 



down to 1.9 magnitude s fainter than the WR s tar in optical R-band observa- 
tions. Soon thereafter, IChurchwell et ahl (jl992l ) made multi-wavelength radio 
observations with the Very Large Array (VLA) and confirmed the nature of the 
emission components, but also did not identify a companion star. A compact 
companion was s uggested to exp l ain th e X-ray and non-thermal radio emission. 
It was not until Williams et al. ( 19971 ) carried out sufficiently high-resolution 



imaging at 2.2 microns with UKIRT that a nearby companion was identified, 3 
magnitudes fainter than the WR star at 2.2 microns and with the luminosity of 
a BO-type star. The companion is 60 mas to the north of the non-thermal radio 
component, in a position consistent with the stellar wind momentum ratio of a 
WN8 and a BO star (see Fig. [1]) and with enough kinetic luminosity to power 
the X-ray emission. 

Subsequently, multi- frequency radio data of WR 147 was used to attempt to 

const rain models of the emission from CWBs (e.g. Skinner et al.|[T99^ : Doughertv et al 



20031 ) . However, because the thermal emission from the stellar wind of the WN8 



star in WR 147 dominates at radio wavelengths, the continuum spectrum of the 
non-thermal emission is not of sufficient signal-to-noise to constrain the models 
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Figure 1. Left: Overlay of a MERLIN 5-GHz image of WR147 (contours) 
and a 2.2/xm UKIRT image (greyscale). The IR image of the WR star was 
ahgned with the peak position of the thermal emission in the southern ra- 
dio component, assumed to be the ionized stellar wind of the WR star. 
The position of the non-thermal emission region (the northern radio com- 
ponent) is consistent with the location of ram pressure balance between the 
winds of the WR sta r and the B-type companion revealed in the IR image 
(IWilHams et al.lll997D . Right: Radio fluxes from WR147 with model fits 
to both the total (solid symbols) and the non-thermal emission component 
(open symb ols) based on 2 D-hydrodynamic models of the WCR and stellar 
winds (from [Pittard et all 12006). Note the large uncertainties in the fluxes, 
particularly at the higher frequencies. 



well. Better model constraints have had to wait for the observations of other 
systems. 

2.1.1 WR 146 - toward better modelling constraints 

At a similar time as the WR 147 observations, iDoughertv et al" ( 19961 ) found ev- 



idence in MERLIN 5-GHz observations for a wide colliding-wind binary system 
in WR 146. This assertion was supported by observations from th e HST that re- 



veale d an 0-type star 168 milli-arcsecond north of the WC star (jNiemela et al.l 



19981 ). Unlike WR147, the radio emission from this system is dominated by 



the non-thermal emission arising from the WCR and so is a more ideal object 
for both very high resolution multi-wavelength imaging and high-quality ra- 
diometry that give more rob ust modelling constrai nts. Initial models have been 
constructed for this system ( O'Connor et al.ll2005l ) (see Fig [5]). 



3 WR 140 - the particle acceleration laboratory 

Both WR 146 and WR 147 are "visual" binary systems. As such, the WCR oc- 
curs far out in the stellar wind of the dense WR star, and the lines-of-sight to the 
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Figure 2. Left: Best-fit models of the radio emission from the WCR in 
WR146, with synchrotron (dotted), thermal (lower dot-dashed) and syn- 
chrotron-|-thermal (dashed). The thermal emission from the WR star wind 
is the higher dot-dashed line. The solid line is the total flux. Right: High- 
resolution images of WR146 from VLA 43-GHz (blue) and VLBA 8.6-GHz 
(red) observations. The radio observations are aligned after accounting for 
proper motion between the two observations. The 43-GHz observations show 
the southern thermal source that is associated with the WR star and the 
northern non-thermal emission from the WCR. The VLBA observation re- 
veals the structure of the WCR. The crosses mark the relative location of 
the two stars, assuming the WR star is coincident with the southern 43-GHz 
thermal source (O'Connor et al., in preparation). 



WCR are typically optically thin. The apparent separation of the stars (asinz, 
where a and i are the orbit size and inclination respectively) in these systems 
suggest that, with an assumed binary mass of ~ 50Mq, the orbital periods are 
a minimum of ~ 350 yrs. In such systems, variations due to the relative motion 
of the stars will all have long timescales - long relative to the time frame of 
typical study! Such systems lack a well-defined "orbit", most particularly the 
separation of the stars. In spite of the apparent success modelling the contin- 
uum spectra and spatial distribution of the radio emission in both WR146 and 
WR147 (see lPittard et al.ll2006l : [O'Connor et al.|[2b05, ). this presents an imped- 
iment to modelling as a number of key system criteria remain ill-defined e.g. 
the UV-radiation density and the wind-kinetic luminosity available within the 
WCR. 

Closer systems, with periods of a few years and orbits that are well con- 
strained through spectroscopic studies, yet are sufficiently wide for optically thin 
lines-of-sight to the WCR through the circumbinary stellar wind plasma are bet- 
ter suited for modelling. In addition to excellent system constraints, those with 
particularly eccentric orbits have variable conditions within the WCR that offer 
an additional probe of the underlying physics. 

The 7.9-yr period, highly eccentric WC7-I-04 binary WR 140 is one such 
system, and is the archetype CWB system. It is well k nown for its dramatic 
variations in emission from near-IR to radio wavelengths ( Williams et al. 1990l : 



Radio observations of massive stars 



5 





0.2 0.4 



Figure 3. Left: 5-GHz VLA radiometry of WR 140 as a function of orbital 
phase for orbit cycles between 1978-1985 (pentagons), 1985-1993 (squares), 
1993-2001 (triangles ), and th e current cycle 2001-2009 (circles ). Op en symbols 
are from the VLA |white fc Beckc ii 119951: [Do ughert y at al.li200ah and solid 
symbols from the WSRT (IWilliams et al.lll99a il994) . Right: A montage of 
8.4-GHz VLBA observations of WR 140 at three different orbital phases that 
demonstrate the rotation o f the WC R as the orbit progresses. The deduced 
orbit is superimposed (Dougherty et al.,,2005.) . 



White k Beckerll 19951). an d also at X -ray energies throughout its orbit (|Zhekov &: Skinneii 



20001 : iPollock et alJ l200"l . l2005l ). The WCR in WR140 experiences significant 
changes in its local environment as the stellar separation varies between ~ 2 AU 
at periastron and ~ 30 AU at apastron. The observed radio emission increases by 
up to two orders of magnitude between periastron and a frequency-dependent 
peak between orbital phases 0.65 and 0.85, followed by a steep decline. This 
behavior repeats almost exactly from one orbit cycle to another, suggesting a 
particle acceleration mechanism that is well-controlled by the orbit of the system 
(Fig. El). 



3.1 The first fully defined CWB orbit geometry 

To fully establish the geometry of the WR140 binary system, particularly the 
orbital inclination (i) and semi-major axis (a), but also the longitude of the as- 
cending nod e (Q), the system mus t be resolved into a "visual" binary. This was 
achieved by iMonnier et al.l (|2004l ) using the Infrared-Optical Telescope Array 
(IOTA) interferometer . Together with other established orbit parameters (e.g. 
Marchenko et al. 20031 ) . families of possible solutions for (i,r2,a) can be deter- 



mined. Until similar observation at other orbital phases are available, VLBI 
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Figure 4. A couple of models of the radio emission from WR 140 at or- 
bital phase 0.837. Note the importance of a high-precision 43-GHz datum to 
distinguishing models. 



observations of the WCR are the only means to determine uniquely i, and hence 
0, and a, from these solution families. 

A 24-epoch campaign of VLB A observations of WR 140 was carried out 
between orbital phase 0.7 and 0.9 ( Doughertv et al.ll2005l ). At each epoch, an 
arc of emission is observed, resembling the bow-shaped morphology expected 
for a WCR. This arc rotates from "pointing" NW to W as the orbit progresses. 
The WCR emission is expected to wrap around the star with the lower wind 
momentum - the O star. In this case, the rotation of the WCR as the orbit 
progresses implies that the O star moves from the SE to close to due E of the WR 
star over the period of the VLBA observations. Thus, the direction of the orbit 
is established. Secondly, the orbital inclination can be derived from the change 
in the orientation of the WCR with orbital phase. Each (i, fi) family provides 
a unique set of position angles for the projected line-of-centres as a function of 
orbital phase. This gives i = 122° zt 5° and f2 = 353° it 3°, leading to a semi- 
major axis of a = 9.0 it 0.5 mas, and the first full definition of orbit parameters 
for any CWB system. From optical spectroscopy a sini = 14.10 it 0.54 AU 
(|Marchenko et al.l l200j ). implying a distance to WR140 of 1.85 ± 0.16 kpc. 



Notably, this estimate is independent of a calibration of stellar parameters e.g. 
absolute magnitude, and implies that the luminosity of the O star is that of a 
super giant. 
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Figure 5. Left: An 8.4-GHz image of Cyg 0B2 #5 from the VLA, that 
clearly shows two components. The bright component is "thermal" and as- 
sociated with the 06-binary system, and the weak component to the NE is 
non-thermal. Right: The radiometry at 4.8 and 8.4 GHz of the bright com- 
ponent, folded with a 6.7-yr orbit period and showing the best-fits at the two 
frequencies of a simple free-free opacity model with a non-thermal source in 
a 6.7-yr eccentric orbit around the binary (Kennedy et al, in preparation). 



3.2 Modelling the radio emission from WR140 

With the geometry and stellar parameters, models of the WCR in WR 140 have 
better observational constraints than in any other known CWB. In addition, 
observations of the thermal X-ray emission that arises from the shocked plasma 
in the WCR provide the best available constraints on the mass-loss rates of the 
two winds as a function of the m omentum ratio of the winds. 

Pittard Doughertvl ( 20061 ) studied models of the radio and X-ray emission 



at orbital phase 0.837, where non-thermal emission dominates the radio spec- 
trum and a good estimate of the characteristics of the relativistic electrons is 
possible. Fig. [3] shows two model fits to the radiometry, with quite different val- 
ues of relative wind-momentum. The models suggest the low-frequency turnover 
at ~ 3 GHz is due to free-free absorption in the winds, since models where the 
Razin effect dominates the low-frequency spectrum place an unacceptably large 
fraction of the shock energy in non-thermal electrons. Development of these 
models continues (see Pittard, these proceedings). 

Observations over a wide radio frequency range, and particularly those 
at the higher and lower frequencies, are key to distinguishing between models 
(see Fig. H|). Currently, radiometry is the only available technique since high- 
resolution radio observations made to date do not have sufficient low brightness 
sensitivity combined with resolution to detect lower brightness emission far out 
in the downstream flow of the WCR that may reveal the asymptotic opening an- 
gle of the WCR, a function of wind-momentum ratio. The current upgrade to the 
Very Large Array (the EVLA) will significantly improve radiometry precision, 
and the broad-banding of VLBI networks e.g. VLBA, cMERLIN will improve 
high-resolution imaging capability, most especially when combined with the data 
from the EVLA. 
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Figure 6. Left: 5-GHz flux of Cyg 0B2 #9, modulated by a 2.35-year 
period, based on a preliminary reduction of archive VLA data. Phase 0.0 was 
arbitrarily set at radio minimum. Right: An 8.4-GHz VLBA observation of 
Cyg 0B2 #9. A bow-shaped WCR is clearly resolved (Dougherty et al., in 
preparation) 



4 Looks like a duck, quacks like a duck - what is it? O-star binaries 

The bulk of observations and modelling of CWB systems has been established 
using radio-bright WR-I-OB binary systems. Given that these are the descen- 
dents of massive O-star binaries, it begs the question of whether the CWB model 
can also explain non-thermal emission amongst O stars. Certainly, of the 16 O- 
star systems that exhibit non-thermal emission (~ 40% of all the radio-d etected 
O stars), at least 70% are established binary systems ( De Becker! l2007l ). Here, 



a couple of systems with good observational evidence for a WCR are given as 
examples. 

In the 6.6-day binary 06-1-06 star syste m Cyg 0B2 # 5, is both a thermal 



and a non-thermal radio component (e.g. iMiralles et al.lll994). similar to those 



observed in WR 146 and WR 147 (Fig. [5]). IContreras et al.l (|l997l ) suggested that 
the non-thermal component is from a WCR between the O-star binary and a 
B-type star ~ 0.9 arcseconds distant. However, most recently Kennedy et al. (in 
preparation) have demonstrated that radio em ission from the putat ive WCR is 
constant and the radio variations first noted bv lMiralles et all ( 19941 ) arise solely 



from the "thermal" component that is associated with the 06 binary. They 
succeed in modelling the variations through orbit-modulated free-free opacity 
along lines of sight to a non-thermal source in a 6.7-yr orbit around the 06-star 
binary, where at minimum the spectral index is closely that of a "classic" stellar 
wind plasma (~ 0.6), and at maximum is flatter (~ 0.25) due to the non -thermal 
contribution to the total emission (Fig. [5]). Recentlv lLkider et al.1 (12009") suggest 
that this model may also explain the X-ray emission. The success of this simple 
model is compelling, though multi-frequency observations may reveal a more 
complex story - as discovered for WR140! 

In Cyg 0B2 ^9, analysis of the radio emission from the system shows 
that the emission is clearly modulated with a 2.35-year period, very similar to 
the behavior of the radio emission from WR 140. VLBA observations of this 
system reveal a bow-shaped non-thermal emitting region, as anticipated for a 
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wind-collision region (Fig. [B]) . Togeth er, this evidence s uggests strongly a WCR 
origin for the emission, and recently i Naze et al. ( 20081 ) succeeded in the long- 
standing challenge of finding spectroscopic evidence of the companion required 
to give rise to the WCR. 

The observational support for the CWB model among 0-type stars is rising 
(see van Loo, these proceedings). The evidence from imaging and/or radiometry 
in these examples of recent radio work is remarkably similar to those of the 
established WR+OB CWB systems - even in the absence of supporting evidence 
for a stellar companion. The evidence of WCRs in non-thermal emitting O stars 
continues t o grow in new obser vations, particularly at higher resolution e.g. 



HD93129A (|Benagha et al.ll2009l 1 



5 Finding new systems 

Identifying stellar companions through the traditional technique of detecting ra- 
dial velocity variations can be extremely challenging. It is strongly dependent 
on the luminosity (and mass) contrast of the two stars in the binary systems, 
and the period and orientation of the orbit. Certainly for many of the currently 
known CWB systems, the periods are long (> lyr) and requ ire long-term d edi- 
cated observing programmes over many years to identify fsee lWiniamil200l for 
examples) . 

Identifying CWB systems is important in a number of leading areas of 
astrophysics. Establishing the binary fraction among massive stars is crucial 
for models of binary star formation, especially in young massive stellar clusters 
where the long-period binary systems may be readily disrupted. Secondly, the 
abili ty of the WC-type WR binaries to produce quite large amounts of dusts 
(e.g. Tuthill et al.lll999l ) leads to the question of whether massive binary systems 



could have provided significant dust production during the era of first stars and 
galaxy formation. 

The CWB systems that we know of today have been identified through at 
least one, and often two of the following: 1) non-thermal radio emission, 2) 
thermal X-ray emission or, 3) variable dust emission, all emission from within 
the WCR. A number of methods that take advantage of these observational at- 
tributes have been successfully used recently and suggest methods for identifying 
binary systems that will not be readily identified using spectroscopic techniques. 
Furthermore, these methods mitigate against the need for long-term monitoring 

to d etect binarie s. 

Clark et"an d^OOS") studied the X-ray properties of the young massive clus- 



ter Westerlund 1 and showed that a number of members have hard {kT > 
2.0 keV) and bright (L.^ ^ 10^^~ ^^ erg s~^) X-ray emission as anticipated from 



a WCR (e.g. IStevens et al.llT99^ ). Interestingly, these stars also exhibit either 



radio continuum s pectra that are too flat to be thermal emission alone (e.g. 



Doughertv &: Clar k 2008 ), and /or they exhi bit an IR excess that has been at 



tributed to hot dust (jCrowther et al. l2006l ). These identifications have con- 



tributed to determinin g a very high bin ary fraction (> 70%) amongst the WR 



stars in Westerlund 1 (Clark et al 



5 n bin a 
200811. 



Perhaps most dramatically, Tuthill et al. ( 20061 ) used diffraction- limited IR 



imaging to readily identify the five 'cocoon' stars in the massive young Quin- 
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tuplet Cluster as WR-binary systems, with each having a characteristic CWB 
dusty "pinwheel" nebula. The difficulty of previous attempts at identifying the 
character of these stars, demonstrates the power of this high-resolution tech- 
nique, certainly for identifying dusty WC-type WR stars. 

Most recently, the ChlcAGO project has been identifying "anomalous" 
ASCA X-ray sources with fluxes in the range 10~^^ — 10~^^ erg s~^ in the 
Galactic Plane. Such sources are not bright enough to be well-known, yet too 
bright to be canonical active stars or AGN. Through follow-up observations with 
CHANDRA (for resolution), 2MASS and spectroscopy (for stellar identification) 
and radio observations a mmibcr of these objects have now been identified as 
colliding- wind systems (G. Anderson &; B. Gaensler, priv. communication). 

6 The future of this research 

High-resolution imaging at radio wavelengths has provided dramatic confirma- 
tion that relativistic electrons are being accelerated in wind-collision regions in 
massive binary systems. The evolution of the WCR revealed in WR 140 shows 
what can be achieved with this technique. However, high quality radiometry is 
still central to models of the radio emission, and other high-energy emission, in 
colliding- wind binaries. 

With the advent of the upgrade to the Very Large Array, the EVLA, it 
will be possible to attain radiometry at least an order of magnitude better than 
currently available. Further, it will be possible to observe the entire continuum 
spectrum from 1 to 50 GHz with unprecedented precision. The Atacama Large 
Millimetre Array (ALMA) is another exciting development that will have the 
ability to image dust emission structures around massive stars, most particularly 
older, cooler dust that exists further from the underlying WC binary systems. 

The ability to readily identify binaries through the radio. X-ray and IR char- 
acteristics of CWBs gives a very powerful tool for extending our understanding 
of the birth-rate and formation environments of massive stellar binaries. This is 
especially timely given the recent discoveries of more massive, optically obscured 
stellar clusters in the Milky Way, especially in the direction of the Galactic Cen- 
tre. 

Acknowledgments. Much of the work described in this brief review was 
done in collaboration with many colleagues. Most particularly, I have had many 
engaging discussions with both Julian Pittard (Leeds University) and Perry 
Williams (IfA, Edinburgh) as part of a number of successful projects related 
to colliding- winds. Also I like to thank a number of very enthusiastic students 
who have worked with me on various CWB projects over the past six years 
- Nick Bolingbroke, Evan O'Connor, Marshall Kennedy and Thomas Johnson. 
Lastly, I'm grateful to Gemma Anderson for sharing some of the results from 
the ChlcAGO project prior to publication. 

References 

Abbott, D. C, Beiging, J. H., ChurchwcU, E., & Torres, A. V. 1986, ApJ, 303, 239 
Abbott, D. C, Bieging, J. H., & Churchwell, E. 1984, ApJ, 280, 671 



Radio observations of massive stars 



11 



Benaglia, P, Dougherty, S. M., Phillips, C, Koribalski, B., Tzioumis, T. 2009, Proc. 

Interferometic view on Hot Stars, to appear in Revista Mexicana de Astronomia 

y Astrofisica Conference Series 
Churchwell, E., Bieging, J. H., van der Hucht, K. A., Williams, P. M., Spoelstra, 

T. A. T., & Abbott, D. C. 1992, ApJ, 393, 329 
Clark, J. S., Muno, M. P., Ncgucrucla, I., Dougherty, S. M., Crowther, P. A., Goodwin, 

S. P., & de Grijs, R. 2008, A&A, 477, 147 
Contreras, M. E., Rodriguez, L. F., Tapia, M., Cardini, D., Emanuele, A., Badiali, M., 

& Persi, P. 1997, ApJ, 488, L153+ 
Crowther, P. A., Hadfield, L. J., Clark, J. S., Negueruela, I., & Vacca, W. D. 2006, 

MNRAS, 372, 1407 
De Becker, M. 2007, A&A Rev., 14, 171 

Dougherty, S. M., Beasley, A. J., Claussen, M. J., Zauderer, B. A., & Bolingbroke, N. J. 

2005, ApJ, 623, 447 

Dougherty, S. M., & Clark, J. S. 2008, in Revista Mexicana de Astronomia y Astrofisica 

Conference Series, Vol. 33, 68 70 
Dougherty, S. M., Pittard, J. M., Kasian, L., Coker, R. P., Williams, P. M., & Lloyd, 

H. M. 2003, A&A, 409, 217 
Dougherty, S. M., & Williams, P. M. 2000, MNRAS, 319, 1005 
Dougherty, S. M., Williams, P. M., & PoUacco, D. L. 2000, MNRAS, 316, 143 
Dougherty, S. M., Wilhams, P. M., van der Hucht, K. A., Bode, M. F., & Davis, R. J. 

1996, MNRAS, 280, 963 
Eichler, D., & Usov, V. 1993, ApJ, 402, 271 

Jardine, M., Allen, H. R., & Pollock, A. M. T. 1996, A&A, 314, 594 

Linder, N., Rauw, G., Manfroid, J., Damerdji, Y., De Becker, M., Eenens, P., Royer, 
P., & Vreux, J.-M. 2009, A&A, 495, 231 

Marchenko, S. V., Moffat, A. F. J., Ballereau, D., Chauville, J., Zorec, J., Hill, G. M., 
Annuk, K., Corral, L. J., Demers, H., Eenens, P. R. J., Panov, K. P., Seggewiss, 
W., Thomson, J. R., & Villar-Sbaffi, A. 2003, ApJ, 596, 1295 

Miralles, M. P., Rodriquez, L. F., Tapia, M., Roth, M., Persi, P., Ferrari- Toniolo, M., 
& Curiel, S. 1994, A&A, 282, 547 

Monnier, J. D., Traub, W. A., Schloerb, F. P., MiUan-Gabet, R., Berger, J.-P., Pedretti, 
E., Carleton, N. P., Kraus, S., Lacasse, M. G., Brewer, M., Ragland, S., Ahearn, 
A., Coldwell, C, Haguenauer, P., Kern, P., Labeye, P., Lagny, L., Malbet, F., 
Malin, D., Maymounkov, P., Morel, S., Papaliolios, C, Perraut, K., Pearlman, 
M., Porro, I. L., Schanen, I., Souccar, K., Torres, G., & Wallace, G. 2004, ApJ, 
602, L57 

Moran, J. P., Davis, R. J., Spencer, R. E., Bode, M. F., & Taylor, A. R. 1989, Nat, 340, 
449 

Naze, Y., De Becker, M., Rauw, G., & Barbieri, C. 2008, A&A, 483, 543 
Niemela, V. S., Shara, M. M., Wallace, D. J., Zurek, D. R., & Moffat, A. F. J. 1998, 
AJ, 115, 2047 

O'Connor, E. P., Dougherty, S. M., Pittard, J. M., & Williams, P. M. 2005, in Massive 
Stars and High-Energy Emission in OB Associations, ed. G. Rauw, Y. Naze, & 
E. R. Blomme, Cosset, 81-84 

Pittard, J. M., & Dougherty, S. M. 2006, MNRAS, 372, 801 

Pittard, J. M., Dougherty, S. M., Coker, R. F., O'Connor, E., & Bolingbroke, N. J. 

2006, A&A, 446, 1001 

Pollock, A. M. T., Corcoran, M. F., & Stevens, I. R. 2002, in Interacting Winds from 
Massive Stars, ed. A. F. J. Moffat & N. St-Louis, Vol. 260, 537—1- 

Pollock, A. M. T., Corcoran, M. F., Stevens, I. R., & Williams, P. M. 2005, ApJ, 629, 
482 

Skinner, S. L., Itoh, M., Nagase, F., & Zhekov, S. A. 1999, ApJ, 524, 394 
Stevens, I. R., Blondin, J. M., & Pollock, A. M. T. 1992, ApJ, 386, 265 



12 



Dougherty 



Tuthill, P., Monnier, J., Tanner, A., Figer, D., Ghez, A., & Danchi, W. 2006, Science, 
313, 935 

Tuthill, P. G., Monnier, J. D., & Danchi, W. C. 1999, Nat, 398, 487 

van der Hucht, K. A., Williams, P. M., Spoclstra, T. A. T., & de Bruyn, A. C. 1992, in 

Nonisotropic and Variable Outflows from Stars, ed. L. Drissen, C. Leitherer, & 

A. Nota, Vol. 22, 249-+ 
White, R. L., & Becker, R. H. 1995, ApJ, 451, 352 

Williams, P. M. 1996, in Radio Emission from the Stars and the Sun, ed. A. R. Taylor 

& J. M. Paredes, Vol. 93, 15-+ 
Williams, P. M. 2002, in Interacting Winds from Massive Stars, ed. A. F. J. Moffat & 

N. St-Louis, Vol. 260, 311~+ 
Williams, P. M., Dougherty, S. M., Davis, R. J., van dcr Hucht, K. A., Bode, M. F., & 

Setia Gunawan, D. Y. A. 1997, MNRAS, 289, 10 
Williams, P. M., van der Hucht, K. A., Pollock, A. M. T., Florkowski, D. R., van der 

Woerd, H., & Wamsteker, W. M. 1990, MNRAS, 243, 662 
Williams, P. M., van dcr Hucht, K. A., & Spoclstra, T. A. T. 1994, A&A, 291, 805 
Williams, P. M., van der Hucht, K. A., & The, P. S. 1987, QJRAS, 28, 248 
Zhekov, S. A., & Skinner, S. L. 2000, ApJ, 538, 808 



